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ABSTRACT 

We carried out extended spectroscopic confirmations of Lya emitters (LAEs) at z — 6.5 and 5.7 
in the Subaru Deep Field. Now, the total number of spectroscopically confirmed LAEs is 45 and 54 
at z = 6.5 and 5.7, respectively, and at least 81% (70%) of our photometric candidates at z = 6.5 
(5.7) have been spectroscopically identified as real LAEs. We made careful measurements of the Lya 
luminosity, both photometrically and spectroscopically, to accurately determine the Lya and rest-UV 
luminosity functions (LFs). The substantially improved evaluation of the Lya LF at z — 6.5 shows 
an apparent deficit from z = 5.7 at least at the bright end, and a possible decline even at the faint 
end, though small uncertainties remain. The rest-UV LFs at z = 6.5 and 5.7 are in good agreement, 
at least at the bright end, in clear contrast to the differences seen in the Lya LF. These results imply 
an increase in the neutral fraction of the intergalactic medium from z = 5.7 to 6.5. The rest- frame 
equivalent width (EWo) distribution at z = 6.5 seems to be systematically smaller than z — 5.7, and 
it shows an extended tail toward larger EWo- The bright end of the rest-UV LF can be reproduced 
from the observed Lya LF and a reasonable EWo— UV luminosity relation. Integrating this rest-UV 
LF provides the first measurement of the contribution of LAEs to the photon budget required for 
reionization. The derived UV LF suggests that the fractional contribution of LAEs to the photon 
budget among Lyman break galaxies significantly increases toward faint magnitudes. Low-luminosity 
LAEs could dominate the ionizing photon budget, though this inference depends strongly on the 
uncertain faint-end slope of the Lya LF. 

Subject headings: cosmology: observations — early Universe — galaxies: formation — galaxies: high- 
redshift 
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Cosmic reionization was a major event in the early his- 
tory of the universe. It is a drastic phase transition of the 
intergalactic medium (IGM) and is closely related to the 
birth of the first galaxies; however, it is still unclear when 
and how reionization occurred. The polarization mea- 
surement of the cosmic microwave background (CMB) 
by Wilkinson Microwave Anisotropy Probe (WMAP) im- 
plies reionization at z = 10.9 ± 1.4 (Komatsu et al. 
2009), and the complete Gunn-Peterson (CP) trough of 
the Sloan Digital Sky Survey (SDSS) QSOs suggests that 
cosmic reionization ended at z ~ 6 (e.g., Fan et al. 2006). 
In addition to these unknowns, it is unclear which ob- 
jects were responsible for the cosmic reionization. Al- 
though QSOs are expected to be the main contributor 
at the bright end of the luminosity function (LF) of ion- 
izing sources, the QSO population alone cannot account 
for all the required ionizing photons (Jiang et al. 2008, 
Willott et al. 2005). Star- forming galaxies such as Ly- 
man break galaxies (LBGs) and Lya emitters (LAEs) 
at the reionization epoch are the only alternatives that 
could dominate the LF at the faint end. 

The LAEs are one class of high- z star- forming galaxies. 
Detecting their strong Lya emissions is feasible even be- 
yond z = 6 using deep narrow band (NB) imaging. Along 
with LBGs, they tell us about early star formation his- 
tory and initial structure formation. In addition, they 
may serve as valuable observational tools for probing the 
cosmic reionization process. Lya photons are sensitive 
to the physical state of the IGM. It is expected that the 
damping wing from the surrounding neutral IGM atten- 
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uates Lya photons so significantly that the Lya emission 
flux will be reduced; therefore, it is naturally expected 
that the Lya LF of LAEs should decline as it traces ear- 
lier times in the reionization epoch (Haiman & Spaans 
1999; Malhotra & Rhoads 2004). The LAE population 
appears to have a similar number density over a long- 
time period from z = 3 to 6 (e.g., Ouchi et al. 2008). 
Thus, if a sharp decline appears in the Lya LF of LAEs, 
it could result from attenuation by the neutral IGM. Con- 
sequently, the observed census of LAEs during the reion- 
ization period could trace the neutral fraction of IGM 
hydrogen, xhi- The census of observable LAEs beyond 
z = 6 is sensitive to xhi of the universe. 

This "Lya test" has an advantage in that it is sensi- 
tive even at Xm > 10~ 3 , which is the upper limit for the 
application of the GP test. In addition, this method uses 
galaxies that are more abundant than QSOs or gamma- 
ray bursts (GRBs), so it will yield a volume-averaged 
estimate for the neutral fraction. Significant variation in 
IGM transmissions among different QSO lines-of-sight 
(Djorgovski et al. 2006) suggests that the reionization 
process is spatially patchy. In the future, it will be pos- 
sible to investigate field-to-field variation in the neutral 
fraction, providing qualitative estimates of this spatially 
patchy reionization process. Alternatively, the change in 
the Lya LF could be caused by galaxy evolution. Distin- 
guishing IGM attenuation from galaxy evolution is cru- 
cial for this test. However, recent systematic surveys for 
LAEs at lower- z revealed that the Lya LF exhibits al- 
most no evolution from z — 3 to z — 6, though the phys- 
ical reason for this is unclear. In any case, this result 
strongly supports the viability of the Lya test, in which 
any decline of such static LAE LF at z > 5.7 would be 
caused by an IGM attenuation. Nevertheless, a part of 
the LF decline should be caused by the intrinsic galaxy 
evolution of LAEs; therefore, it is important to compare 
the rest-UV LFs based on the LAE sample at the same 
time as the Lya test is made, because the rest-UV con- 
tinuum flux is not sensitive to the neutral IGM. 

Several model predictions of the Lya LF at the reion- 
ization epoch agree that the amplitude of the LF de- 
creases as xhi increases (McQuinn et al. 2007; Haiman 
& Cen 2005; Le Delliou et al. 2005; Dijkstra et al. 2007b; 
Mesinger & Furlanetto 2007; Kobayashi et al. 2007; Iliev 
et al. 2008; Dayal et al. 2009). The H n bubbles of 
bright LAEs clustered in the overdense regions would 
effectively overlap, creating a larger H n bubble with a 
high ionization fraction, which would significantly reduce 
the Lya flux attenuation. As a result, the bright LAEs 
that are expected to reside in a high-density environment 
should be readily observed, whereas faint LAEs are more 
severely attenuated. This results in a smooth decrease in 
the amplitude of the LF toward higher cchi- To compare 
with these model predictions, we have to observe a large 
sample of high-z LAEs. 

In our previous papers (Kashikawa et al. 2006: K06, 
Shimasaku et al. 2006: S06), we have presented spectro- 
scopic confirmations of 17 and 34 LAEs at z — 6.5 and 
5.7, respectively, in the Subaru Deep Field (SDF); see 
also Kodaira et al. (2003) and Taniguchi et al. (2005) 
for LAEs at z = 6.5. The sample consisted of objects 
with excess flux in narrow- band NB921 (A c = 9196 A, 
FWHM=132 A) and NB816 (A c = 8150 A, FWHM=120 



A) images. We found that the Lya LF at z = 6.5 reveals 
an apparent deficit at least at the bright end, correspond- 
ing to ~ 0.75 magnitudes fainter, compared with that 
observed at z — 5.7. The decline in the Lya LF from 
z = 5.7 to 6.5 could imply a substantial transition in 
the IGM ionizing state between these epochs, suggesting 
that reionization was not complete at z = 6.5. The trend 
was also confirmed in the LAE sample at z w 7 (lye et 
al. 2006). The Lya luminosity density of LAEs did not 
change from z — 3 to z — 5.7, and it gradually decreased 
from z = 5.7 to 7 (Ota et al. 2008). The abrupt drop in 
luminosity density from z = 5.7 to 7.0 did not seem to 
be caused by continuous galaxy evolution. This gradual 
decline recalls the IGM attenuation, though we have few 
spectroscopically identified galaxies at z = 7. 

Turning to the faint end of the LF below L(Lya)= 
5 x 10 42 , the amplitude difference between our 
photometrically- and spcctroscopically-dctcrmincd esti- 
mates was too large to constrain the faint end. Our faint 
spectroscopic sample at both z — 6.5 and 5.7 was still so 
small that we could not reliably identify a significant dif- 
ference in the LF between these two epochs at the faint 
end. We neither determined the faint-end slope of the 
LF nor constrained the true contribution of LAEs to the 
entire photon budget required for full reionization. 

In this paper, we describe our extended spectroscopic 
confirmations of LAEs at z = 6.5 and 5.7 after K06 and 
S06 to determine more accurately the faint end of the 
LAE LFs close to the reionization epoch. Spectroscopic 
confirmations enhance the reliability of the LAE sample, 
eliminating possible contaminations from emission-line 
galaxies at lower z, and improving the derivation of the 
Lya luminosity with precise redshifts. A large number of 
spectroscopically identified LAEs enable us to construct 
more accurate LFs in both the Lya luminosity, which is 
sensitive to the neutral IGM, and the rest-UV luminos- 
ity, which is insensitive to it. The larger spectroscopic 
sample also provides better statistical estimates of the 
contamination rate, increasing the reliability of our pho- 
tometric LAE sample. Purely photometric detection of 
LAE candidates might have non-negligible ambiguities 
in estimating both Lya luminosities, which is sensitive 
to the neutral IGM, and the rest-UV luminosity, which 
is insensitive to it, especially at high-z beyond z = 6.5, 
where it is difficult to detect them even with broadband 
imaging (Hibon et al. 2010; Tilvi et al. 2010; Ota et al. 
2010). An increase in the spectroscopic sample size could 
also improve the stacked spectrum, which may provide 
insights into either the internal dynamics of the LAE or 
the IGM properties (K06). 

The Subaru Suprime-Cam has a very large field of view 
(FOV; 27 x 34 arcmin 2 ). Wide-field imaging increases 
the chance of discovering rare objects, such as the most 
distant galaxies (lye et al. 2006). Furthermore, it is gen- 
erally predicted in almost all models of the reionization 
process that an overlapping H II bubble at the end of 
the reionization could be as large as 10 physical Mpc 
(Barkana & Loeb 2004; Wyithe & Loeb 2004; Furlan- 
etto et al. 2004), which corresponds to the FOV size of 
the Subaru Suprime-Cam. Thus, wide FOV observations 
are essentially required to achieve an unbiased picture of 
the universe at the reionization period. Our LAE sam- 
ples were obtained from a general blank field without 
resorting to amplification of gravitational lensing by fore- 
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TABLE 1 

Summary of Spectroscopic Observations 



Observational run Instrument date(UT) Tj ntog (ksec) seeing(arcsec) Nmagk 21 N65 b Ns7 c 



2006 d FOCAS Apr. 26 7.2—9.0 0.5—1.0 2 2 

2007 FOCAS Mayl9— 22,24 10.8—16.2 0.5—1.2 5 15 4 

2008 DEIMOS Apr.30, Mayl 5.4—10.8 0.7—1.1 4 12 

2009 DEIMOS Apr.26— 27 10.8—16.2 0.7—0.8 3 8 10 
2009 e DEIMOS Apr. 24— 25 10.8—12.6 0.7—0.8 2 2 1 
2010 f FOCAS Mar. 18— 21 18.0 2 3 
Total 28 20 



a The total number of masks. 

b The total number of identified LAEs at 2 = 6.5. 
c The total number of identified LAEs at 2 = 5.7. 

d These were identified in the course of other observational program in Nagao et al. (2007). 

°Thcsc were identified in the course of other observational program in L. Jiang ct al. (2011, in preparation). 

f These were identified in the course of other observational program in M. Iye et al. (in preparation). 



ground clusters, providing reliable statistics about their 
LFs. High-z surveys rendered by gravitational lensing, 
which are accessible to low-luminosity sources (Richard 
ct al. 2008, Stark et al. 2007, Hu et al. 2002), are com- 
plementary to our survey. 

This paper is organized as follows: In Section 2, we 
review our photometric LAE sample and initial spectro- 
scopic identifications, presented in our previous studies. 
In Section 3, we describe our new spectroscopically iden- 
tified LAE sample at z = 6.5 and 5.7. In Section 4, wc 
accurately derive the Lya and rest-UV luminosities of 
our LAE sample, comparing the Lya and UV LFs be- 
tween z — 6.5 and 5.7. Comparisons of the rest-frame 
equivalent width (EWo) and stacked Lya emission pro- 
file are presented in Section 5 and Section 6. In Section 
7, we present a new method for deriving the rest-UV LF 
of LAEs, using the Lya LF and EWo distribution, to 
constrain the photon budget required for reionization. A 
summary of the paper is provided in Section 8, with some 
discussion of the implications for reionization based on 
our results. 

Throughout the paper, we assume cosmology param- 
eters: O m = 0.3, Ha = 0.7, and H — 70 h 70 km s _1 
Mpc -1 . These parameters are consistent with recent 
CMB constraints (Komatsu et al. 2009). Magnitudes are 
given in the AB system. 

2. SUMMARY OF THE PREVIOUS LAE SAMPLE AT 
Z = 6.5 AND Z = 5.7 AS OF 2006 

The sample selection of photometric LAE candidates 
at z = 6.5 and z — 5.7 was presented in Taniguchi et al. 
(2005) and S06, respectively. The details of observation, 
photometry and color selection were presented in those 
papers; here, we briefly discuss about our photometric 
selection of the LAE sample. The sample was based on 
flux-excess objects in narrow-band NB921 and iV.B816 
images, compared with very deep broadband images of 
the SDF (Kashikawa et al. 2004). We selected LAE 
candidates at z = 6.5 with definite NB excesses down 
to the 5cr limiting magnitude of NB921 = 26.0, i.e., 
z' - NB921 > 1.0 and z' - NB921 > (z! - NB921) 3a , 
where (z' — NB921) 3a is the 3a error in the color of 
z 1 — NB921, and the very red color in broad-bands, i.e., 
(i' - z' > 1.3 and z < i' 2a - 1.3) or (z > i' 2a - 1.3), 
where i' 2(7 are defined as 2a limiting magnitudes of i', 



to reduce contamination from foreground emission-line 
galaxies. We also applied the no detection (ND; < 3a) 
criteria in all bands blueward of the dropout band, i.e., 
B, V, and R. In Taniguchi et al. (2005), we found 58 pho- 
tometric candidates at z = 6.5 down to iV£?921 = 26.0 in 
the effective survey region of 876 arcmin 2 ; spectroscopic 
confirmations of 17 candidates were presented in K06. 
In S06, we selected LAE candidates at z = 5.7 with NB 
excesses with i' - NB816 > 1.5 down to AB816 = 26.0 
and very red color in broad-bands, i.e., (R— z' > 1.5 and 
R < R2&) or (R > i?2<r), where i?2<r are defined as 2a lim- 
iting magnitudes of R. We also applied the ND (< 2a) 
criterion in the B and V bands. The 89 photometric can- 
didates were found at z — 5.7 down to NB816 = 26.0 
in the same survey region; 34 of them were confirmed by 
spectroscopy. We estimated the detection completeness, 
defined as the ratio of detected LAE candidates to all 
the LAEs actually present in the universe. This should 
be corrected in evaluating the LF as a function of NB 
magnitudes by counting artificial objects distributed on 
the real NB images. The detection completeness was 
evaluated as > 0.8 at NB816 < 25.0 and - 0.75 at 
NB816 = 26.0, and as > 0.75 at NB921 < 25.0 and 
- 0.45 at NB921 = 26.0, for the z = 5.7 and 6.5 sample, 
respectively. We quantitatively distinguished LAEs from 
nearby emitters based on their asymmetric emission-line 
profile, using the "weighted skewness" indicator (K06). 
The comoving survey volume was as large as 2.17 x 10 5 
Mpc 3 and 1.80 x 10 5 for z = 6.5 and z = 5.7, respectively. 
It should be noted that the LAE samples at these two 
epochs were extracted from the same field using the same 
photometric procedure, similar survey volume, and sim- 
ilar selection criteria. We carefully determined the NB- 
excess criteria to provide almost the same EW thresholds 
(EWo > 10A) for both LAE samples. 

3. NEW SPECTROSCOPIC IDENTIFICATIONS FROM LAE 

SAMPLE 

3.1. Spectroscopic Observations 

We carried out further spectroscopic observations with 
the Subaru FOCAS (Kashikawa et al. 2002) and the Keck 
II DEIMOS (Faber et al. 2003) over the last four years 
(2006 — 2010). This spectroscopic campaign is summa- 
rized in Tabic 1. 

The FOCAS observations were made with a 300-linc 
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TABLE 2 

Spectroscopic Properties of z = 6.5 LAEs 
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a The object IDs are those of T05, except ID— 60 and 61, that are not listed in the photometric catalog of T05. 

b The redshift was derived from the wavelength of the flux peak in an observed spectrum assuming the rest wavelength of Lya to 
be 1215A. These measurements could be overestimated in the case of significant damping wings by IGM. Also, the observed peak 
position was slightly shifted redward due to instrumental resolution. See Figure 13. 

c The observed line flux corresponds to the total amount of the flux within the line profile. The slit-loss was corrected. 

^No dust absorption correction was applied. 

c Corrected for instrumental broadening. 

f Observed with FOCAS (F) or DEIMOS (D). 

g Thc coordinates of T05 ID-47 were incorrectly appeared in T05. 

mm -1 grating and an 058 order-cut filter. The spectra 
covered 5400-10, 000 A, with a pixel resolution of 1.34 A. 
The 0". 6- wide slit gave a spectroscopic resolution of 7.1 
A (R ~ 1300). The spatial resolution was 0".3 pixel -1 
with 3-pixel on-chip binning. We also allocated slits for 
A^i?921-strong (z' — NB921 > 1) emitters, irrespective 
of their (i' — z') color, as a LAE criterion to test our 
selection criteria. Some extra slits were allocated for our 
candidates in another observation by Nagao et al. (2007), 
in which they used a 175-line mm -1 Echcllc grating, an 
SDSS z' filter as an order-cut filter, and 0".83 slit, giving 
R ~ 1500. Extra slits were also allocated for observations 
by M.Iye et al. (in preparation), in which they used the 
VPH900 grism, 058 order-cut filter, and 0".8 slit, giving 
R ~ 1500. We obtained spectra of standard stars Hz 44, 
Feige 34, and BD+28 for flux calibration. The data were 
reduced in a standard manner. 

For the DEIMOS observations, we used an 830-line 
mm -1 grating and a GG495 order-cut filter. The central 
wavelength was set to A c = 8100 A. The slit width was 
1".0 with 0.47 A pixcr 1 , giving a resolving power of ~ 

3600. The wavelength coverage was ~ 5000 10, 000 

A, depending on the position in the mask. We also allo- 
cated slits for strong iVi3921 emitters. Some extra slits 
were allocated for our candidates in another observation 



led by L. Jiang et al. (2011, in preparation), in which 
they used almost the same instrumental setup, except 
that A c = 9239 A. We obtained spectra of standard 
stars BD +28 4211 and Hz 44 for flux calibration. The 
data were reduced with the spec2d pipeline 15 for DEEP2 
DEIMOS data reduction. During the 2010 run, the flex- 
ure compensation system (FCS) of DEIMOS was broken. 
Unexpected flexure at this time was adequately corrected 
by shifting the spectral images based on the instrumental 
rotation angle, which is sensitive to the flexure change, 
of each exposure. 

All spectroscopic observations were taken at high 
enough resolving power to distinguish single Lya emis- 
sion from [O Ii] doublets (rest-frame separation of 2. 78 A); 
however, in practice, it is hard to discriminate between 
the features, given very faint emissions. To quantita- 
tively distinguish LAEs from nearby emitters, we used 
a line asymmetry estimator, called weighted skewness, 
S w , as in K06. Skewness is a popular statistic defined 
as the third moment of the distribution function. In 
K06, we found that this estimator is sensitive to asymme- 

15 The data reduction pipeline was developed at University of 
California, Berkeley, with support from National Science Founda- 
tion grant AST 00-71048. 
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TABLE 3 

Spectroscopic Properties of z = 5.7 LAEs 



ID a 


NAME 


z b 


/ s P ec (Lya) c L s P ec (Lya) d FWHM e 




OBS. f 








(10- 18 ergs s- 1 cm- 2 ) (10 42 /i TO 2 ergs s" 1 ) (A) 


(A) 





34549 


SDF 


J132418. 1+271639 


5.702 


9.56 


3.38 


16.9 


9.97 + 1.99 


D 


36849 


SDF 


J132435. 8+271709 


5.738 


6.92 


2.49 


12.0 


32.8 ± 11.7 


D 


45159 


SDF 


J132520.2+271842 


5.701 


10.5 


3.72 


11.1 


6.09 ± 3.28 


D 


61394 


SDF 


J132456.5- 


-272212 


5.761 


39.6 


14.4 


7.09 


3.45 ± 0.47 


F 


62511 


SDF 


J132419.6- 


-272228 


5.665 


15.0 


5.23 


4.81 


6.52 ± 2.07 


D 


65899 


SDF 


J132438.0- 


-272315 


5.687 


15.0 


5.29 


5.34 


13.8 + 2.13 


D 


66162 


SDF 


J132528.8- 


-272317 


5.688 


16.4 


5.78 


8.01 


12.0 + 2.56 


D 


70773 


SDF 


J132347.1- 


-272414 


5.696 


6.77 


2.39 


11.0 


3.35 ± 2.86 


D 


71751 


SDF 


J132344.8- 


-272427 


5.720 


10.4 


3.70 


9.96 


5.30 + 1.27 


F 


75550 


SDF 


J132447.7- 


-272512 


5.708 


14.9 


5.28 


6.63 


4.74 ± 3.44 


D 


81382 


SDF 


J132451.4- 


-272626 


5.697 


9.05 


3.20 


8.55 


6.20 ± 2.00 


D 


83092 


SDF 


J132441.3- 


-272649 


5.747 


15.6 


5.62 


5.50 


5.70 + 5.22 


F 


89624 


SDF 


J132424.8- 


^272812 


5.714 


13.1 


4.67 


8.36 


4.08 ± 5.53 


F 


91179 


SDF 


J132417.8- 


-272833 


5.712 


13.3 


4.74 


14.0 


8.31 ± 1.35 


D 


98040 


SDF 


J132432.9- 


-273009 


5.674 


25.9 


9.06 


11.0 


7.62 ± 1.80 


F 


98461 


SDF 


J132428.3+273012 


5.697 


11.3 


3.99 


3.99 


3.65 ± 1.81 


F 


120018 


SDF 


J132340.6+273447 


5.684 


38.8 


13.6 


11.3 


10.0 + 2.00 


D 


124905 


SDF 


J132352.9+273549 


5.694 


6.88 


2.43 


5.07 


6.93 ± 4.93 


D 


138608 


SDF 


J132450.9+273839 


5.669 


19.7 


6.88 


11.4 


4.56 ± 10.3 


F 


138624 


SDF 


J132357.4+273836 


5.677 


8.30 


2.60 


11.0 


5.04 + 2.75 


D 



''ID in the NB816-detected catalog of S06 

3 The redshift was derived from the wavclc 
be 1215A. These measurements could be overestimated in the case of significant damping wings by IGM. Also, the observed peak 
position was slightly shifted rcdward due to instrumental resolution. Sec Figure 13. 

c The observed line flux corresponds to the total amount of the flux within the line profile. The slit-loss was corrected. 

d No dust absorption correction was applied. 

e Corrected for instrumental broadening. 

f Observed with FOCAS (F) or DEIMOS (D). 



try and carefully determined the critical value to distin- 
guish LAEs from nearby emitters based on our A5921- 
and iV.B8i6-excess sample. The S w values of foreground 
emitters never exceeded S w = 3, which we set as the 
critical S w value to distinguish LAEs from foreground 
emitters. We should note that the threshold 

Sw — 3 ap- 
plied in this study is only valid for LAEs at z = 5.7 and 
6.5. We have no guarantee that the same critical value 
can be used for LAEs at different z, where, in princi- 
ple, different intergalactic attenuation could change the 
degree of asymmetry in the Lya emission. We serendipi- 
tously identified some LAEs (three for z = 6.5 and eight 
for z = 5.7) from strong NB921 and NB816 emitters that 
did not meet our color selection criteria. We found that 
most of them had very close neighbors in the images, 
which prevented accurate aperture photometry. In sum- 
mary, we identified 28 and 20 additional LAEs at z = 6.5 
and 5.7, respectively 16 . The spectra of these additional 
LAEs are shown in Figures 1 and 2, and their spectro- 
scopic properties are summarized in Tables 2 and 3 for 
z = 6.5 and 5.7, respectively. 

We did not detect Ny AI240, which is the only ac- 
cessible strong high-ionization metal line indicative of 
AGN activity, from any of these newly identified LAEs. 
We have some candidates (two for z — 6.5 and six for 
z = 5.7) that exhibited only a faint, poorly fit "single" 
line emission signal, so we could not verify its asymmetry, 

16 SDF J132417.9+271746 and SDF J132521. 1+272712, classi- 
fied as single-line emitters in Taniguchi et al. (2005), and SDF 
J132417.8+272833 (ID 91179), classified as a single-line emitter in 
S06, were found to be LAEs in the deep spectroscopic observations 
performed in this study. 



which is the only key diagnostic for distinguishing Lya 
emission from other emissions. We did not include these 
single-line emitters in the LAE samples. We had some 
additional candidates (7 for z = 6.5 and 10 for z = 5.7) 
for which no emission signal was detected, even using 
8 — 10m class telescopes. These might be intrinsically 
too faint, or some could be transient objects. Table 4 
summarizes our current spectroscopic identifications of 
LAEs at z — 6.5 and 5.7 in the SDF and Figure 3 shows 
the NB-magnitude histogram of spectroscopic identifica- 
tions. At this time, 90% and 74% of the photometric 
candidates at z = 6.5 and 5.7, respectively, have been 
followed by spectroscopy. Uncertain "single" and "no- 
detection (ND)" candidates are of course only dominant 
at the very faint end of the sample, so we consider this 
to be our current observational feasibility limit. De- 
spite this uncertainty, most of the photometric candi- 
dates have been spcctroscopically identified, and only six 
candidates remain without follow-up spectroscopy in the 
z = 6.5 sample. 

3.2. Spectroscopic Sample 

Now the total numbers of spectroscopically confirmed 
LAEs in the SDF are 45 and 54 at z = 6.5 and 5.7, respec- 
tively. Figure 4 shows the line-peak wavelength distribu- 
tion of spectroscopically confirmed LAEs. The distribu- 
tions of line peaks show an apparent blueward deviation 
against the NB transmission curve, which was already 
confirmed in our previous studies. This is due to the fact 
that the asymmetric LAE profile, with a broad red wing 
and sharp blue cutoff and a discontinuous Lyman contin- 
uum break, contributes more to the NB flux when a line 
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Fig. 1. — Spectra of 28 spectroscopically confirmed LAEs at 
2 = 6.5. ID number and redshift are shown in the upper left corner 
of each panel. Spectra taken by DEIMOS were smoothed with a 
3-pixel boxcar. The vertical dotted line indicates the center of the 
Lya emission line. The vertical scale is marked in 0.1 X 10 — 18 erg 
s _1 cm -2 A -1 , and a scaling factor to obtain the correct scale 
appears in the right corner in each panel. The sky spectrum is 
overplotted on the left bottom panel with an arbitrary flux scale. 

peak lies at the blue side of the transmission curve. The 
evaluations of sample completeness and contamination 
rates are listed in Table 4. At least 81% (70%) of our 
photometric candidates at z — 6.5 (5.7) have been spec- 
troscopically identified as real LAEs. Because "single" 
and "ND" candidates depend on uncertain assumptions, 
we evaluated two extreme values, i.e., we assumed that 
all of these uncertain candidates were really LAEs, and 
then assumed that none of the candidates was really a 
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Fig. 2. — Same as Figure 1, but for 20 LAEs at z = 5.7. 

LAE. By taking the average of these two cases, our pho- 
tometric sample completeness factor, determined as the 
ratio of the number of true LAEs to the number of objects 
that meet our selection criteria (= (1 — (CT))/(CP), see 
Table 4), was evaluated to be (1 - 0.105)/0.935 = 0.957 
and (1 - 0.180)/0.860 = 0.953 for z = 6.5 and z = 5.7, 
respectively. These spectroscopic results indicate that 
our photometric sample has high reliability, indicating 
our color selection process works well. [O in] and [O n] 
emitters identified by spectroscopy based on their dou- 
blet signatures, were removed from the photometric sam- 
ple, whereas LAEs that were serendipitously found by 
spectroscopy but not listed in the original photometric 
sample were included in the following analysis. 

4. LYq AND REST-UV LUMINOSITY FUNCTIONS 

4.1. Deriving Lya Luminosity 

We derived Lya luminosity both photometrically and 
spectroscopically. The latter would be more accurate 
than the former, but hardly gives an accurate measure- 
ment of faint continuum flux in most cases, though these 
are eventually found to be fairly consistent with one an- 
other. We measured luminosities more carefully in this 
than the previous study (K06) to accurately determine 
Lya and rest-UV LFs. We applied the following proce- 
dure, not only to the newly identified LAE sample but 
also to the previously identified objects by recalculating 
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TABLE 4 

Summary of Spectroscopic Identifications 



z 


N p hot a 


Lya 


Ha 


[O in] 


[O ii] 


single 13 


ND C 


wo/spec d 


Lya S orcnd C 


Lya tot ai f 


CP min g 


CP h 

max 


C- I max 


CT min j 


6.5 


58 


42 





1 





2 


7 


6 


3 


45 


93% 


94% 


19% 


2% 


5.7 


89 


46 








4 


6 


10 


23 


8 


54 


85% 


87% 


30% 


6% 



a Photomctric candidates 

^These possess only a single poor-quality symmetric emission feature with small S w < 3, in which we cannot determine whether they are Lya emission 
or [O n] emission (K06). 
c No detection of any emission features 

^Candidates without spectroscopic followup observations 

c LAEs screndipitously discovered in spectroscopic observation. They are not included in the photometric candidates. 
^Total number of LAEs spectroscopically identified — Lya+Lya sercnc i 

g Samplc completeness rate provided that all singles arc foreground objects — Lya/Lyatotal 
^Sample completeness rate provided that all singles arc LAEs — (Lya+singlc)/(Lya tota i4-singlc) 

"Sample contamination rate provided that all singles and NDs are foreground objects — ([O lll],[0 ll]+singlc+ND)/(N p hot-wo/spcc) 
J Samplc contamination rate provided that all singles and NDs arc LAEs — [O lll],[0 ll]/(N p h t-wo/spcc) 
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Fig. 3. — Magnitude histogram of spectroscopically identified 
LAEs. 

the Lya and rest-UV luminosities. 

First, for the photometric estimate of Lya line flux, the 
Lya line (/line! erg s _1 cm -2 ) and rest-UV continuum 
(/ c ; erg s -1 cm~ 2 Hz -1 ) fluxes were evaluated from NB 
and BB magnitudes {ttinb and tubb) as follows: 



m NB ,BB + 48.6 = -2.51og 



Jo Lya (fc + flinc)T N B : BBdv/y_ 



IT, 



NB.BB 



dv/v 



fi) 



where v^ ya is the observed frequency of Lya, and Tnb 
and Tb b are the transmission bandpasses of the NB and 
BB filters, respectively. The (NB, BB) filter combina- 
tions in the above formula are (7VB921, z') and (7VB816, 
z') for z = 6.5 and 5.7, respectively. We assumed that 
the spectral energy distribution (SED) of LAEs had a 
constant / c and <5-function Lya emission profile, whose 
line- width was much smaller than VLya- The estimate of 
the Lya line flux was not largely affected, even when tak- 
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Fig. 4. — Observed wavelength distribution of the Lya line peak 
for our spectroscopic sample at z = 5.7 (left) and 6.5(right). The 
transmission curves of ./V.B816 (left) JV.B921 (right) are overplottcd 
as dashed lines. 

ing account of a finite line width as of observed < 30A, 
which made ~ 16% Lya flux change at most at the red- 
dest wavelength of the NB-band. The observed value 
was adopted for VLya when spectroscopically identified; 
otherwise, the central frequency, v c , of the NB filters was 
used. We used the actual filter transmission curves con- 
volved by telescope, instrument and atmospheric trans- 
missions for Tnb.bb- The importance of such a care- 
ful analysis of filters' non-square transmission curve has 
been pointed out in Gronwall et al. (2007). If an object 
was not detected in the BB filter, ttibb was replaced by 
the la limiting magnitude of BB, and f c was forcibly set 
to when f c < 0. 

In spectroscopic measurements of Lya line flux, we cor- 
rected for slit loss. The target was regarded as almost 
perfectly centered in the slit because the typical position- 
ing error of the slit on the target object was as small as 
0".l rms for both FOCAS and DEIMOS spectrographs. 
Thus, we assumed that the slit flux loss depended only 
on the seeing size, target size, and slit width. We cal- 
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Fig. 5. — Comparison of the Lya fluxes measured by spectra, 
fapec, with those inferred from photometry, fphot i for our spectro- 
scopic LAE sample at z = 6.5. The solid line represents a one- 
to-one correspondence between f spe c and fphot- The errors were 
estimated based on the sky rms fluctuation of each spectrum for 
(spec and errors in magnitudes for fphot- 

culatcd a possible slit flux loss based on the seeing size 
and slit width for each spectroscopic observation. An in- 
trinsic Lya radial light profile of each object was simply 
assumed to be a Gaussian profile with FWHM measured 
on the NB image deconvolved by the seeing size (0".98) 
of the image. The slit lengths were always so much longer 
(> 8".0) than the object size that slit loss along the spa- 
tial direction was negligible. The resultant silt loss was 
corrected for the spectroscopically measured Lya flux for 
each object. The average slit loss was evaluated to be 
24.1 ± 10.1% and 30.3 ± 12.5% for z = 6.5 and 5.7, re- 
spectively. The difference in slit flux loss between these 
two samples was minor. 

Comparisons of the Lya fluxes measured photometri- 
cally and spectroscopically for the spectroscopically iden- 
tified LAE sample are shown in Figures 5 and 6. The 
correspondence is good except in a few cases. Note that 
the photometric Lya flux may have non-negligible un- 
certainty when rriBB is only given by an upper limit, as 
appears often in the z = 6.5 sample. We hereafter use 
spectroscopically measured Lya fluxes for the spectro- 
scopically identified LAEs, and photometrically inferred 
Lya fluxes for the remaining uncertain (referred to as 
"single", "ND", and "wo/spec" in Table 4) objects. 

4.2. Lya Luminosity Function 

We estimated the acceptable Lya LF range specified 
by the upper and lower limits. The upper limit was esti- 
mated assuming that all the uncertain ( "single" , "ND" , 
and "wo/spec" in Table 4) photometric candidates, for 
which photometrically evaluated Lya luminosities were 
adopted, were really LAEs. The lower limit was esti- 
mated assuming that all the uncertain candidates were 
not LAEs, i.e., using only the pure spectroscopically 
identified LAE sample. To derive the Lya LF, the detec- 
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Fig. 6. — Same as Figure 5, but for LAEs at z = 5.7. 

tion completeness, which decreases at fainter NB magni- 
tudes, was corrected in the same way as outlined in K06. 
In both the upper and lower limit estimates, we corrected 
for this detection completeness by number weighting ac- 
cording to the NB magnitude. 

Figure 7 provides a comparison of the cumulative Lya 
LFs between z = 6.5 and 5.7. The red-shaded and blue- 
shaded regions indicate acceptable LF ranges for z = 6.5 
and 5.7, respectively. Compared with our previous es- 
timate in K06, the acceptable LF ranges are now more 
sharply determined. The given uncertainties are by the 
Poisson errors, shown for some average data points be- 
tween upper limit and lower limit data points, indicat- 
ing that Poisson errors are dominant at the bright end, 
whereas the spectroscopic uncertainties are dominant at 
the faint end. Taking into account the corrections for 
completeness factor (95.7% for z = 6.5 and 95.3% for 
z = 5.7) evaluated in the previous section, both LFs 
could increase by a factor of ~ 1.05, although this is 
smaller than the Poisson errors. Note that we derived the 
apparent Lya luminosity uncorrected for either dust ex- 
tinction or self-absorption, which is evident on the blue- 
side cutoff of the emission line. 

We fitted a Schechter function, <j>{L)dL = (j)*(L/L*) a 
exp(—L/L*)dL/L*, to the average data points between 
the upper limit and lower limit data points. The y 2 was 
minimized with a single grid search in the two param- 
eters, L* and </>*, for fixed slopes of a = —1.7, —1.5, 
and —1.3. The larger value of the Poisson errors or the 
difference between the upper limit and the lower limit 
was used for the error of each data point when fitting a 
Schechter function. The derived best-fit parameters are 
listed in Table 5. 

After spectroscopic confirmations of a large amount of 
our LAE samples at both z = 6.5 and 5.7, we found 
that the evaluation of the Lya LF at z = 6.5 was appar- 
ently deficient compared with that of z — 5.7, at least at 
the bright end, where the entire samples at both epochs 
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Fig. 7. — Comparison of the cumulative Lya LFs of LAEs at 
z = 6.5 (red-shaded region) and z = 5.7 (blue-shaded region). We 
estimated the acceptable Lya LF ranges as specified by the upper 
and lower limits. The upper limit was estimated assuming that 
all the uncertain (referred to as "single", "ND", and "wo/spec" in 
Table 4) photometric candidates are really LAEs, and the lower 
limit was estimated assuming that all the uncertain candidates are 
not LAEs, i.e., using only the pure spectroscopically identified LAE 
sample. In both the upper and lower limit estimates, we corrected 
for the detection completeness by number weighting according to 
the NB magnitude. Error bars evaluated by the Poisson errors are 
shown in some average data points between the upper limits and 
lower limits. The short-dashed lines (red for z = 6.5 and blue for 
z = 5.7) show the fitted Schechter LFs in the case of a = —1.5. 
As a comparison, the green long-dashed line shows the Lya LF at 
z = 6.5 from Ouchi et al. (2010), and the green dot-dashed line 
shows that of Hu et al. (2010). 



TABLE 5 

Best-fit Schechter Parameters for LAE LF at 
z = 6.5 and 5.7 OF the SDF 



Sample 


a 


L* 


<f>* 




(fix) 


log(/hf 2 ergs s" 1 ) 


log(//if Mpc" 3 ) 


z = 6.5 


-1.7 
-1.5 
-1.3 


42.8218;!° 
42.7618;!° 

42.7018:18 


-3 40+ ' 22 
o os +0.20 
o - zo -0.20 
, r, n +0.20 
,r> - zu -0.18 


z = 5.7 


-1.7 
-1.5 
-1.3 


43 i2+ u - u « 
43 n9+ 06 
42 Q4+°- oe 


o 7z ,+U.UB 

-3 56+ 08 

, 04+0.10 



have been spectroscopically identified. A possible decline 
in the z ~ 6.5 LF appeared even at the faint end. The 
L* difference between z = 6.5 and z = 5.7 decreased 
slightly to ~ 0.65 magnitudes from the previous esti- 
mate of ~ 0.75 magnitudes, given a — —1.5. A major 
cause of this change could be the significant increase in 
the number of spectroscopically confirmed LAEs and the 
more careful measurements of luminosity performed in 
the present study. 

4.2.1. Uncertainties in the LF estimates 



Two uncertainties appeared in the photometric esti- 
mate of the Lya luminosity: the photometric error and 
the redshift ambiguity. In the previous section, we as- 
sumed that the peak of the Lya emission line was the 
central wavelength, A c , of the NB filters for the photo- 
metric sample. We carried out a Monte Carlo simula- 
tion to see how the associated uncertainties would affect 
the resultant LF. In the simulations, Gaussian random 
photometric errors both in BB and NB were assigned 
to the measured magnitudes of the photometric sample, 
and redshifts were assigned so that their redshift dis- 
tributions matched the observations, as shown in Fig- 
ure 4. The la photometric error was estimated from 
the background fluctuation and the flux Poisson noise. 
The process of recalculating the LF and deriving the 
best-fit Schechter parameters was repeated many times. 
We found rms fluctuations of cr(log(L*)) = 0.041 and 
<r(log(</>*)) = 0.066 for z = 6.5 and cr(log(L*)) = 0.023 
and a(\og((/>*)) = 0.027 for z = 5.7 for a = -1.5, sug- 
gesting that these uncertainties have only a small impact 
on the result. The photometric sample selection prior to 
the luminosity measurements might be affected by errors 
in the photometric catalog. We also carried out a Monte 
Carlo simulation to see how the photometric error in the 
catalog would affect the color selection and resultant LF, 
by assigning random errors to the measured magnitudes 
of detected object catalogs in all bands. We did not use 
spectroscopic estimate of the Lya flux in the simulation, 
and estimated the Lya LFs photometrically. The rms 
fluctuations were found to be cr(log(£*)) = 0.104 and 
cr(logO*)) = 0.119 for z = 6.5 and a(log(L*)) = 0.098 
and <r(log(</>*)) = 0.094 for z = 5.7 for a = -1.5. We 
note that more than 70% of our LAE samples have been 
spectroscopically identified, whereas we have to add the 
redshift ambiguity in all the sample in this simulation; 
therefore, the uncertainties of the LF estimates are ap- 
parently overestimated. 

Given an observed NB flux, the photometric Lya flux 
tended to be overestimated, and the continuum flux 
tended to be underestimated when the actual line peak 
wavelength was smaller than the A c due to a strong Ly- 
man break in the continuum, as well as to the asymmet- 
ric Lya line profile. This trend would be expected if the 
transmission curve of the NB filter was nearly a top-hat 
shape; however, we used NB filters with almost Gaussian- 
shaped transmission curves, as shown in Figure 4. For 
these filters, photometrically estimated line flux would 
maintain the observed NB excess by decreasing when the 
line peak shifted farther away, whether it was redder or 
bluer, from the A c . As a result, the photometric line 
flux was largely underestimated when an emission line 
was really located in the redder part of the NB trans- 
mission, whereas in the blue portion, an underestimate 
due to the low transmission was balanced by an overesti- 
mate due to the continuum Lyman break. We confirmed 
these trends by numerical experiment. As the observed 
line-peak distribution deviated to blue, the systematic 
error in the photometric estimate of Lya line would be 
small. In our Monte Carlo simulation, we found that 
the systematic error caused by this trend was as small 
as a(log(L*)) = 0.04 and 0.06 for z = 5.7, and 6.5, re- 
spectively. We also performed a Monte Carlo simulation 
to investigate any possible distortion that the discrep- 
ancy, as seen in Figure 5 and 6, between spectroscopically 
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measured and photometrically inferred Lya luminosities 
could have on the result. When Gaussian random error 
with the same scatter as in Figure 5 and 6 was assigned to 
each Lya luminosity, the best-ht Schechter parameters, 
given a = —1.5, only changed by <r(log(L*)) = 0.023 and 
cr(log(>*)) = 0.042 for z = 6.5 and er(log(X*)) = 0.042 
and <r(log((/>*)) = 0.070 for z = 5.7, which were negligi- 
ble. 

4.2.2. Comparisons with other studies 

The Schechter parameters derived in this study are al- 
most identical to previous estimates by S06 for z = 5.7, 
whereas at z = 6.5, the values in this study fell be- 
tween the previous estimates of the upper and lower lim- 
its. These parameters are also consistent, within the er- 
rors, with independent studies by Ouchi et al. (2010) 
and Ouchi et al. (2008) for z = 6.5 and 5.7. The es- 
timates at z = 6.5 given by Hu et al. (2010) differed 
by the largest amount from our estimate, especially at 
faint luminosities. Their estimated Lya LF was almost 
a factor of three (five) lower than estimates of K06, 
Ouchi et al. (2008) and Ouchi et al. (2010) at their 
faintest bin at logL(Lya)^ 42.8 at z = 6.5 (5.7). Hu 
et al. (2010) claimed that the difference might be mainly 
caused by a large contamination in the photometric sam- 
ple, if present; however, the present study, based on a 
large number of spectroscopic confirmations, completely 
disallows this interpretation. As shown in Section 3, the 
contamination of our sample was as low as < 20%, eval- 
uated by the spectroscopic results, and the sample com- 
pleteness was as high as > 85% for both the z = 6.5 and 
5.7 samples. The photometric LAE selection of Ouchi 
et al. (2008) and Ouchi et al. (2010) nearly matches 
this study, so their sample is probably also less affected 
by contamination. The reason for the difference in the 
two evaluations at the faint end is unclear because the 
data reduction, photometry, and the LAE selection differ 
slightly. However, we can suggest one plausible explana- 
tion. The Lya LF at z = 6.5 by Hu et al. (2010), which 
is basically based on their spectroscopic sample, is very 
close to that of our spectroscopic estimate (K06). Spec- 
troscopic confirmation at the faint end is generally so dif- 
ficult that our previous spectroscopic sample apparently 
lacked completeness at the faint end. The present study, 
based on spectroscopic observations with deeper typical 
integration times (10 ksec, see Table 1) than Hu et al. 
(2010), overcame the problem and revealed that most of 
these faint unidentified objects were real LAEs. As al- 
most half of the faint photometric sample of Hu et al. 
(2010) were not spectroscopically identified at z = 5.7, 
the completeness difference at the faint end could explain 
the significant difference in resultant Lya LFs. Nonethe- 
less, it is interesting to note that both studies, more or 
less, detected a difference in the Lya LFs between z = 6.5 
and 5.7 

4.2.3. Significance of the LF difference and the cosmic 

variance 

To illustrate the significance of the LF difference be- 
tween z = 6.5 and 5.7, we plotted the error contours for 
our Schechter-parameter fits in Figure 8. The confidence 
levels of the fitting were computed based on the larger 
of the Poissonian error or the difference between the up- 
per limit and the lower limit. Figure 8 reveals that the 
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Fig. 8. — Error ellipses of the best-fit Schechter parameters <f>* 
and L* of Lya LF given a fixed a = —1.5. The lower ellipse (red) 
is for LAEs at z = 6.5, and the upper ellipse (blue) is for z = 5.7. 
The inner and outer solid ellipses are the 1 a and 3cr confidence 
levels, respectively. 

(L* , (j>*) error ellipses at fixed a = —1.5 for z = 5.7 and 
6.5 do not overlap each other; that is, the difference in 
LF between z = 5.7 and 6.5 is significant at almost the 
3(7 level. This is also the case for any a. The differ- 
ence in L* is more significant than that in c/>* . Based on 
Somerville et al. (2004), we evaluated the cosmic vari- 
ance of our LAE samples. We assumed a one-to-one cor- 
respondence between LAEs and dark haloes, and used 
their predictions at z = 6. With our comoving survey 
volume of 2.17 x 10 5 h^ Mpc 3 and a number density 
of 2.76 x 10~ 4 (1.94 x 10" 4 ) /if Mpc" 3 for the upper 
(lower) limit estimate, we obtained a cosmic variance of 
~ 32%. We also estimated a variance of ~ 20% for the 
z = 5.7 LAE sample. As shown by the error bars in 
Figure 8, the 3er error circles for the two epochs overlap 
each other when cosmic variance is included; however, 
our upper limit estimate still differed from the z = 5.7 
result at the 2a level. Although the difference in LF 
was large at the bright end and small at the faint end, 
we forcibly attempted to fit the LF at z — 6.5 for the 
Schechter function at z — 5.7 with either fixed <jj* or 
L*, given a = -1.5. We obtained log L* = 42.90±g:?| 
given a fixed (j>* , corresponding to a 24% decrease in 
the Lya luminosity, L*(z = 6.5) = 0.76L*(z = 5.7), 
and we obtained log (</>*) — —3.74^0 24 given a fixed 
L*, corresponding to a 34% decrease in LF amplitude, 
4>*(z — 6.5) = O.660*(z = 5.7). These values are compa- 
rable to the cosmic variance, indicating that the observed 
LF difference might be caused by the cosmic variance. 
We cannot rule out this possibility, but we discuss it fur- 
ther in Section 8. 

4.3. Rest-UV Continuum Luminosity Function 

We derived the rest-UV continuum LFs of our LAE 
sample at z = 6.5 and 5.7. In Section 4.1, the flux of 
the rest-UV continuum (/ c ) was simultaneously derived 
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once the Lya line flux (/ii nc ) was determined by either 
spectroscopy or photometry. The effective wavelength 
of the derived UV luminosity was 1250 - 1270A for 
both z = 5.7 and 6.5. The detection completeness in the 
z' band, which corresponds to the rest-UV flux, should 
be corrected when calculating the UV LF; however, our 
LAE samples were basically selected in NB magnitude. 
It is impossible to evaluate the detection completeness 
in the z' band for the NB-selected sample; therefore, we 
should note that it is inevitable that the derived rest-UV 
continuum LF may be affected, especially at the faint 
end of the LF, by the difference in completeness of the 
NB and z' bands. As in our previous study, the cor- 
rection was made based on the detection completeness 
in the NB filter evaluated in the previous section. To 
overcome the problem, we take an alternative approach 
to derive the rest-UV LF of LAEs in Section 7. No cor- 
rection has been applied for dust. LAEs are generally 
recognized to be young, less massive galaxies with lit- 
tle dust (E(B - V) < 0.05 at z ~ 6 (Ono et al. 2010) 
and even at z ^3 (Gronwall et al. 2007)), though old 
and massive stellar populations are found in some LAEs 
(Finkelstein et al. 2009). 

Figure 9 shows the rest-UV continuum LF of our LAE 
sample. As in the Lya LF, the red- and blue-shaded 
regions mark the acceptable UV LF ranges for z = 6.5 
and 5.7, respectively, and Poisson errors are shown for 
some data points. The vertical dotted lines indicate the 
corresponding limiting magnitudes in the z 1 band. Our 
LF measurements at magnitudes fainter than Muv — 
—20.24 (3 a) may be uncertain because the correspond- 
ing z'-band magnitudes are no longer reliable. The de- 
rived rest-UV LFs presented in this study are nearly con- 
sistent with our previous studies discussed in K06 and 
S06, though the amplitudes are slightly lower than be- 
fore. We confirmed the result of our previous study, that 
the rest UV LF does not change between z = 6.5 and 
5.7, at least at the bright end of M uv < -20.5. We 
fit the Schechter function to these data points down to 
Mjjv = —20.24, beyond which our measurements were 
largely affected by incompleteness of source detection. 
Because data points were limited to the bright end, we 
used a fixed slope of a = —1.5. The best fit parameters 
were M* v = -21.720^;^ and log(^) = -4.15^ 
for z = 6.5, and M* v = -21.845±?;^ and log(0*) = 

— 4.30^q'7q for z = 5.7. Figure 10 shows the error con- 
tours, which overlap at the la errors. This is in clear con- 
trast to the difference seen in the Lya LF. It should be 
noted that the rest-UV continuum flux is not sensitive to 
the neutral IGM. This result strongly supports the inter- 
pretation that the difference in Lya LF between z = 6.5 
and 5.7 is caused by the IGM attenuation. However, the 
best-fit parameters of rest-UV LFs include large uncer- 
tainties, and the large error contours in Figure 10 suggest 
that our LAE sample is still insufficient to strongly con- 
strain the parameters. Possible distortions in the rest- 
UV LF caused by photometric error were also evaluated 
by the Monte Carlo simulation as in Section 4.2. As- 
suming a = —1.5, the la fluctuations were found to be 
a(M* v ) = 0.253 and <r(log(<£*)) = 0.136 for z = 6.5 and 
a{M uv ) = 0.154 and a (log (<£*)) = 0.087 for z = 5.7. 
These fluctuations are smaller than the fitting errors of 
the best-fit Schechter parameters, though larger than 
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Fig. 9. — Comparison of the rest-UV continuum cumulative 
LFs of LAEs at z = 6.5 (red-shaded region) and z = 5.7 (blue- 
shaded region). The upper and lower limits of the shaded regions 
are determined in the same way as in the Lya LF estimate. The 
vertical lines indicate the limiting magnitudes in the z' band at 
M uv = -19.05, -19.80 and -20.24 for la, 2a, and 3a given 
EWo = at z = 6.5, respectively. The rest-UV LF measurements 
at magnitudes fainter than the 3c limiting magnitude may be un- 
certain due to incompleteness in the z' band; we show the rest-UV 
continuum LFs down to the la limiting magnitude just for ref- 
erence. We did not use data points at > Mjjy = —20.24 when 
fitting the Schechter function. Error bars are calculated from Pois- 
son errors. As a comparison, the rest-UV LF of the LAE sample 
at z = 5.7 (magenta short-dashed line) and z = 3.1 (magenta long- 
dashed line) evaluated by Ouchi et al. (2008), the range of the 
rest-UV LF of the LBG sample at z ~ 6 (green shaded region) 
determined by various studies (see Figure 11 of Bouwens et al. 
2007), and the rest-UV LF of z ~ 3 LBG (Reddy et al. 2008; green 
long-dashed line) are shown. The black dot-dashed line shows the 
estimate of the UV LF at z = 5.7 using Lyo LF with a = —1.5, 
assuming a universal EWo-UV luminosity relation (see Section 7). 

those of Lya LFs. 

In Figure 9, we overplotted another rest-UV contin- 
uum LF estimate (with a simple extrapolation at the 
faint end) of the LAE sample at z — 5.7 by Ouchi et al. 
(2008), which was almost consistent with our measure- 
ment at z = 5.7 at Muv < —20.5. The green-shaded re- 
gion is the range of the UV LF of LBGs sampled at z ~ 6 
determined by various studies, based on Bouwens et al. 
(2007), Figure 11. The coincidence of UV LFs of LBG 
and LAE at z ~ 6 was confirmed by S06 and Ouchi et 
al. 2008. This indicates that most LBGs at z ~ 6 display 
Lya emissions, though the UV LF measurements remain 
uncertain for both in LAE and LBG. It is quite inter- 
esting to note that the UV LFs of these two populations 
are consistent at z ~ 6 because they generally seem to 
show different evolutionary trends at z < 6. In Figure 9, 
we also present the UV LF of LAEs at z = 3.1 (Ouchi et 
al. (2008); magenta long-dashed line) and that of LBGs 
at z ~ 3 (Reddy et al. (2008): green long-dashed line). 
The UV LF of the LAEs seems to show an increase in L* 
at higher z, whereas that of the LBGs seems to have the 
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Fig. 10. — Error ellipses of the best-fit Schechter parameters <f>* 
and L* of the rest-UV continuum LFs of LAEs at z = 6.5 (red 
upper ellipse) and z = 5.7 (blue lower ellipse), assuming a fixed 
a = —1.5. The inner and outer solid ellipses are the 1<t and 3<r 
confidence levels, respectively. The confidence levels of the fit were 
computed based on the larger of the Poisson error or the difference 
between the upper limit and the lower limit. 

reverse trend, showing a decrease in L* . A high fraction 
of LAEs among LBGs in the early universe is naturally 
expected if a LBG appears as a LAE during its initial 
starburst phase, when it is still dust-free (Shapley et al. 
2001). The scenario is supported by measurements of the 
LAE dark halo mass of ~ 10 11 M Q , which is system- 
atically smaller than that of LBGs (Ouchi et al. 2010) 
at all epochs. Shimasaku et al. (2006) and Ouchi et al. 
(2008) reached almost the same conclusion concerning a 
high fraction of LAEs among LBGs at z ~ 6, whereas 
Dow-Hygelund et al. (2007) claimed that the fraction 
of LAEs among LBGs at z ~ 6 was ~ 30%, which was 
almost the same as at z ~ 3. Stark et al. (2011) es- 
timated that the LAE fraction of luminous LBGs with 
-21.75 < M uv < -20.25 at z ~ 6 was as small as 20%, 
though their sample was restricted to strong Lya emis- 
sions with EWo > 25A. A much higher fraction of LAEs 
should be expected at smaller EWo from their EW distri- 
bution. Henry et al. (2010) suggested another possibility 
based on their blind multislit spectroscopic search for 
LAEs: that the z-dropout color selection might miss a 
certain fraction of LAEs that have blue (i — z) colors due 
to strong Lya emission. Otherwise, the UV LF of LBGs 
might be underestimated due to dust extinction, which is 
generally small for LAEs. We cannot make any conclu- 
sive argument about the LAE fraction of LBGs and its 
evolution solely from our observed rest-UV LF; however, 
one should be careful when comparing among different 
observations because the fraction depends on the thresh- 
old EW applied when selecting the LAEs, as well as on 
the rest-UV luminosity. 

5. EW DISTRIBUTION 



A more straightforward observable probe of reioniza- 
tion is the equivalent-width distribution of LAEs, which 
can be derived from the Lya and UV continuum fluxes. 
These are sensitive and insensitive to the neutral IGM, 
respectively. It was very hard to measure the rest-frame 
equivalent width (EWo) m spectroscopic data because 
most LAEs were too faint to accurately measure their 
continuum flux directly from spectra. Instead, the EWo 
was calculated using narrow- and broad-band photom- 
etry. EWo was photometrically estimated in some pre- 
vious work, though most of the EWo values were only 
lower limits because the continuum emission was not de- 
tected in broad-band images. In contrast, most LAEs in 
the present study were actually detected in the z'-band 
image by virtue of the extreme depth of the SDF images. 
The EWo were reliably determined in those cases. 

Figure 11 compares the EWo distribution of the LAE 
samples for z — 6.5 and 5.7. The detection completeness 
was corrected by number weighting according to the NB 
magnitude, as outlined in Section 4.2. Here, we did not 
correct for the absorption of the blue side of the Lya 
emission due primarily to the interstellar medium (ISM) 
absorption inside the galaxy. K06 showed that the blue- 
side line profile of Lya emission of the composite spec- 
trum of LAEs at z — 6.5 was simply explained by spec- 
tral broadening, which meant that the blue side of the 
observed Lya emission was almost completely absorbed. 
Therefore, if the ISM absorption is corrected, the evalu- 
ated EWo would almost double. For several LAEs that 
remained undetected in the z' band (< la), we used the 
lower limit of EWo by replacing the z'-band magnitude 
with the lcr limiting magnitude of z' = 27.79. The sys- 
tematic errors caused by the process will be discussed 
later. 

In Figure 11, we also plotted the EWo distribution 
of lower-z LAEs at z = 3.1, 3.7, and 5.7 extracted from 
Ouchi et al. (2008). The EW distribution of our z = 5.7 
LAE sample was almost consistent with other low-z sam- 
ples, whereas the EWo at z = 6.5 seemed to be systemat- 
ically smaller than that at z = 5.7. This appeared more 
clearly in the cumulative distribution of EWo shown in 
the bottom panel of Figure 11. We evaluated a possible 
distortion in the EWo distribution caused by photomet- 
ric error and the redshift ambiguity by the Monte Carlo 
simulation, as discussed in Section 4.2. The la fluctua- 
tions of the EWo distribution of the simulation, shown 
as the shaded regions in the bottom panel of Figure 11, 
ensured that our results were not seriously affected by 
these uncertainties. S06 suggested that an estimate of 
the EWo distribution at z = 5.7 based on i' and NB816 
magnitudes was largely different from one based on z' 
and NB816 magnitudes. In the former estimate, the Lya 
line enters both the i' and NB816 bands, as it does at 
z = 6.5, which is based on z' and NB921 magnitudes. 
The blue dashed line in the bottom panel of Figure 11 
indicates the EWo distribution at z — 5.7 when EWo is 
calculated from i' and NB816 magnitudes. It approaches 
the EWo distribution at z — 6.5, with an extended tail 
at large EWo- This might be caused by the fact that 
the deeper limiting magnitude of the i' band imposed a 
stronger constraint on the continuum flux than did the 
z'-band, which only provided a lower limit on EWo when 
undetected; therefore, deeper BB magnitudes generally 
tended to yield higher EWq. Excluding objects that were 
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Fig. 11. — Top: the differential EWo fraction distribution of 
the LAE sample at Z = 6.5 (red line) and z = 5.7 (blue line). 
Comparisons with low-z LAEs at z = 3.1, 3.7, and 5.7 from Ouchi 
et al. (2008) are also shown. Bottom: the cumulative EWo fraction 
distribution of the LAE sample at z = 6.5 (red line) and z = 5.7 
(blue line). The orange and cyan shaded regions are shown as 
la fluctuations in the EWo distribution caused by photometric 
errors and redshift uncertainties. The blue dashed line indicates the 
cumulative EWo fraction distribution of the LAE sample at z = 5.7 
given that EWo was calculated from i' and NB816 magnitudes. 



Fig. 12.— UV luminosity-EWo relation at z = 5.7 (left) 
and 6.5(right). Filled circles present the spectroscopic sample, 
and open circles show the photometric sample. The objects with 
EWo >240 A are plotted at the EWo =240 A. Error bars denote 
uncertainties caused by photometric errors and redshift ambigu- 
ity, except for data points that have large uncertainties due to 
no detection in the BB band (only the upper limit of Mjjy and 
the lower limit of EWo are shown as arrows). The dotted curves 
represent the EWo at fixed Lya luminosity; from top to bottom, 
5,2,1,0.5,0.2,0.1 xlO 43 ergs" 1 . 



not detected in the z'-band actually reduced the number 
of objects with high EWo, but did not make a substan- 
tial change in the overall EWo distributions shown in 
Figure 11, i.e., EWo at z — 6.5 seems to be system- 
atically smaller than at z = 5.7. However, there is no 
way to estimate the EWo distribution at z = 6.5 in the 
case of using a BB band, in which the Lya line does not 
enter, because we do not have any sufficiently deep BB 
bands at longer wavelength than the z'-band. We can- 
not completely rule out the possibility that the difference 
seen in the EWo distribution between z = 5.7 and 6.5 is 
attributed to the difference whether the Lya line enters 
into the BB band or not. A very deep J-band photom- 
etry to estimate EWo at z — 6.5 is required to make a 
more fair comparison of the EWo distribution. 

The median EWo at z — 6.5 is 74 A , which is smaller 
by 15A than that at z = 5.7. This trend may be easily 
understood from the resultant LFs of these two epochs. 
The UV continuum LFs are almost identical, whereas 
Lya LF at z = 6.5 is deficient compared with that at 
z — 5.7. The trend was also indirectly suggested from 
a possible FWHM difference in lines seen by Hu et al. 
(2010), for which it was difficult to derive EWo due to 
the shallow photometric data. Fontana et al. (2010) sug- 
gested that the absence of prominent La lines in their 
LBG sample at z ~ 7 compared with a conservative EW 
distribution at lower-z. 

Another interesting feature of this plot is that the EWo 
distribution at z = 6.5 has a remarkable extended tail 
toward larger EWo compared with z — 5.7. Such an 
effect was also found in the LBG sample (Stark et al. 
2007). One explanation for this might be contamina- 



tion by Population Ill-dominated galaxies. It is gener- 
ally thought that the universe was first metal-enriched 
by this first generation of stars. As these Population III 
stars were born in extreme metal-free conditions, they 
presumably had a top-heavy initial mass function. Given 
their exceptionally high effective temperatures, they are 
expected to produce a very hard spectrum, enough to 
ionize the primordial He gas. The main characteristics 
of the predicted SED are the presence of a large-EW 
Lya emission line, due to the strong ionizing flux, and 
strong He 1 and He 11 recombination lines, due to spec- 
tral hardness (Schaerer 2002). Therefore, Population Ill- 
dominated galaxies are expected to appear among high-z 
LAE samples with large EWs. Several LAEs at z = 6.5 
have EWo larger than ~ 300 A, which cannot be at- 
tained by the usual Population II synthesis. These ob- 
jects are plausible Population III candidates because of 
their extraordinarily large EW at high z, though they 
are expected to dominate at higher z (7 < z < 15), ac- 
cording to most model predictions (Johnson et al. 2008, 
Yoshida et al. 2007). One alternative possible origin for 
the large Lya EW is a contribution from AGNs. Other- 
wise, the large Lya EW might be the result of scattering 
in a clumpy, dusty interstellar medium (Neufeld 1991; 
Hansen & Oh 2006; Finkelstein et al. 2009). Follow-up 
NIR spectroscopy to detect the Hell emission signal is the 
only promising way to confirm Population Ill-dominated 
galaxies. However, it should be noted that such an ex- 
tended tail toward larger EWo can, more or less, be pro- 
duced artificially by the uncertainty in BB flux. 

Figure 12 shows the EWo-UV luminosity relation. For 
faint z'-band magnitudes (< 27.79; ler), both EWo and 
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Myy are only provided as lower limits, indicated as ar- 
rows in the figure. Note that the z 1 bandpass directly 
corresponds to the UV continuum luminosity at z = 5.7, 
whereas the Lya flux also contributes to the z'-band flux 
at z = 6.5. A clear vertical lower-limit sequence of Mjjv 
can be seen at z = 5.7 in Figure 12, and not at z = 6.5. 
As reported in previous lower- z studies (S06; Stark et al. 
2010; Vanzella et al. 2009; Ouchi et al. 2008; Ando et 
al. 2006), we saw an apparent deficit of LAEs with large 
EWo at bright UV magnitudes, and the maximum EWo 
increased with lower UV luminosity, which is expected if 
low- luminosity galaxies are less obscured by dust. This 
study confirmed the trend at z = 5.7, and almost the 
same tendency was found at z — 6.5, as well. 

6. LYq PROFILE OF THE COMPOSITE SPECTRUM 

The Lya emission-line profile is also a useful reion- 
ization signature, in principle, because the neutral IGM 
imposes a damping absorption feature on it (Dayal et al. 
2008, Dijkstra et al. 2007a). A difference in Lya line 
profiles between two epochs around the reionization pe- 
riod, might suggest changing IGM opacity. Though there 
are many model predictions of Lya profiles during the 
reionization epoch, numerous factors such as star forma- 
tion rate (SFR), internal kinematics, inflow/outflow, and 
source clustering easily affect the profiles. On the obser- 
vational side, individual spectra are usually too noisy to 
verify the line profile, so several spectra must be stacked 
to obtain good average line-profile features. As in our 
previous study (K06), the composite spectrum was made 
using the procedure outlined below. We now have 45 
and 54 LAE spectra at z = 6.5 and 5.7, respectively, 
though at different spectroscopic resolutions. First, we 
removed the spectra with the poorest instrument resolu- 
tion. Then, each spectrum was smoothed with a Gaus- 
sian kernel chosen to produce a common instrument reso- 
lution of FWHM~ 6.4 A, which was practically measured 
from the FWHM of sky lines near the Lya emission for 
z = 5.7 and 6.5, respectively. Each spectrum was shifted 
so that the line peak wavelength, which is the only fea- 
ture used to estimate the redshift, was at the rest wave- 
length of 1215.67A. The spectra were rebinned to a com- 
mon pixel scale, and then coadded by taking the average 
with scaling and weighting based on their line flux, using 
a 3<T clipping to eliminate sky-subtraction residuals. 

Figure 13 compares the final composite spectrum be- 
tween z = 6.5 and 5.7. Both composite spectra revealed 
an apparently asymmetric profile with an extended red 
wing. We found no significant differences between these 
two composite spectra, even when comparing the com- 
posite spectra made only from Lya-bright or Lya-faint 
objects. This conclusion is consistent with Ouchi et al. 
(2010) and Hu et al. (2010). The composite spectrum 
at z = 6.5 seems to have a slight excess over z — 5.7 at 
the red wing tail from 1217A to 1220A, though it is much 
smaller than the la error of the difference in the compos- 
ite spectra. This could be caused by a slight broadening 
in the Lya emission- line profile itself from z — 5.7 to 
6.5, as suggested by Ouchi et al. (2010). Otherwise, this 
could be caused by scaling to systematically smaller Lya 
flux at z = 6.5, compared to that of z = 5.7 (as seen on 
the EWo distribution), which increases the UV contin- 
uum level at z = 6.5. 

The similarity in Lya emission profiles might indicate a 
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Fig. 13. — Comparison of the composite spectrum between z = 
6.5 (red) and 5.7 (blue). The lower panel shows the residual of 
subtraction of z = 5.7 from 6.5. The dotted line shows the la 
error of the difference of the composite spectra. 

lack of significant IGM opacity evolution from z = 5.7 to 
6.5. However, using the composite spectrum of Lya emis- 
sion to constrain the reionization has a critical problem 
of its own. As mentioned above, the wavelength of the 
Lya line peak has to be used to measure the redshift, and 
is often systematically offset from the rest frame. Shap- 
ley et al. (2003) found the kinematic offset implied by 
the relative redshifts of Lya emission and low-ionization 
interstellar absorption lines in an LBG sample at z ~ 3. 
The offsets vary significantly as a function of Lya emis- 
sion strength from 800 km s _1 to 480 km s — 1 , though 
these shifts could be smaller at high z, where the emit- 
ting halos are substantially less massive, and might not 
power such strong winds as at z ~ 3. This can be ex- 
plained in terms of the properties of large-scale outflows 
(Adelberger et al. 2003, Mas-Hesse et al. 2003, Westra 
et al. 2005). Recently, McLinden et al. (2011) were the 
first to detect the [O III] emission line from two LAEs 
at z ~ 3, and they also found apparent velocity off- 
sets. Such systematic offsets in Lya emission lines from 
the systemic redshift would increase the uncertainty in 
measuring redshifts using only the Lya peak. Conse- 
quently, when making a composite spectrum using sev- 
eral LAE spectra, these systematic offsets would dilute 
the line profile, critically preventing us from revealing an 
accurate line profile. Also, if the field-to-field variances 
of the neutral fraction and the size of H II bubble were 
large even at the same redshift, the composite line profile 
would be diluted. It seems hard to derive a conclusive 
constraint on reionization from line profiles alone. The 
most promising method for overcoming the problem is 
to measure the systemic redshift of LAEs by detecting 
nebular emissions (Adelberger et al. 2005, McLinden et 
al. 2011), which is still technically difficult at z > 5. This 
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wind effect could allow the Lya emission line to emerge 
at wavelengths where the GP optical depth is reduced, 
transmitting the Lya flux directory to the observer even 
at the reionization epoch (Dijkstra et al. 2010). Mc- 
Quinn et al. (2007) found that a 400 km s _1 redshift in 
the Lya line did not have a large effect on their conclu- 
sions regarding the effect of reionization on the Lya LF 
and clustering at xhi < 0.4, and they concluded that this 
wind effect would not seriously affect the "Lya test" in 
the late phase of reionization. 

We would also like to note that the Lya radiative trans- 
fer is complicated by the geometry and kinematics of the 
ISM and IGM. The relative geometries of interstellar H 
I and H II regions significantly affect resonant scatter- 
ing, which can either suppress or enhance the Lya line 
(Chariot & Fall 1993, Neufeld 1991). The resonant scat- 
tering in the IGM induces a change in the frequency of 
Lya photons, which can also cause a mass-dependent 
redward shift of the Lya line peak, even in the absence 
of a galactic wind (Zheng et al. 2010). Dust attenua- 
tion, if any, also significantly reduces the Lya emission. 
A stacking analysis would, more or less, overlook these 
variations among galaxies. 

In K06, we showed that our composite Lya line pro- 
file at z = 6.5 could be realized by both the reioniza- 
tion model, in which we included the attenuation of a 
GP damping wing from outside the H n bubble, and 
the galactic wind model, which has another broadly ex- 
tended Gaussian component in the line profile. Although 
additional spectroscopic data have improved the quality 
of the composite spectrum, the spectral resolution of our 
composite spectrum is still too low to distinguish between 
these two models. Higher resolving power provided by 
larger telescopes will be required to constrain the model 
more strongly. 

7. CONTRIBUTION OF LAES TO THE REIONIZATION 
PHOTON BUDGET 

The integration of the observed UV LF at the faint 
end provides an estimate of the luminosity density and, 
thus, of the photon budget of reionization. Therefore, an 
accurate UV LF estimate based on the LAE sample will 
constrain the LAE contribution to the photon budget. 
Unfortunately, the UV LF measurements at magnitudes 
fainter than Mjjy = —20.24 (3a) may be critically un- 
certain because of incompleteness, which cannot be cor- 
rected due to the difference in the completeness of the 
NB and z'-band data. Here, we take another approach 
to derive the UV LF of LAEs using both the Lya LF and 
the EWo-UV luminosity relation, as shown in Figure 12. 

The EW distribution function is well approximated 
by an exponential function (Gronwall et al. 2007; Stark 
et al. 2011), though its e-folding width, w, might de- 
pend on UV luminosity. Several previous works provide 
derivations of the EWo-UV luminosity relation of LAE 
at lower-z (e.g., Stark et al. 2010; Vanzella et al. 2009; 
Ouchi et al. 2008; Ando et al. 2006). Their conclusions 
are roughly consistent with an apparent deficit of large 
EWo at bright UV magnitudes and a maximum EWo 
that increases with lower UV luminosity. The relation 
does not show a clear evolution across the 3 < z < 6 
redshift range (Ouchi et al. 2008, Cassata et al. 2011), 
and our result at z = 6.5 also roughly follows this trend 
(Figure 12), though we derived our result from a rel- 



atively small sample. Among these previous studies, 
Stark et al. (2010) derived a reliable EWo-UV luminos- 
ity relation based on a large spectroscopic LAE sample 
at 3 < z < 7. Based on this relation, the characteris- 
tic e-folding width, w, was empirically determined to in- 
crease with fainter UV magnitudes: w = &QMjjv + 1440, 
which gives w = 60A at Mjjy = —23 and w — 300A at 
Mjjv = —19, respectively. This equation almost traces 
the maximum EWo values as a function of Mjjv ■ Figure 
12 of Stark et al. (2010) shows an apparent deficiency in 
LAEs with low EWo and low UV luminosity, though it 
could be caused by their spectroscopic detection limit. 
Hence, we assumed that the peak of the exponential dis- 
tribution was always at EWo = 0, independent of Mjjv- 
A cutoff EW, below which LAEs were not included in 
the sample, was set to 10A , based on the LAE selec- 
tion criterion m z — m^B > 1, and Equation (1). With 
this EWo probability distribution, we conducted Monte 
Carlo simulations to derive Muv, given a Lya luminos- 
ity and EWo- The Lya LF measurements included an 
accurate correction for completeness; therefore, the UV 
LF was estimated free of incompleteness even at the faint 
end when assuming the EWo distribution, as described 
above. 

The dot-dashed line in Figure 9 shows the estimate 
of the UV LF at z = 5.7 with this method. Here, we 
assumed that the faint end of the Lya LF at z = 5.7 
was a = —1.5. Note that we did not use the UV LF 
at z = 6.5 to constrain the photon budget because the 
Lya LF at z = 6.5 might be affected by neutral IGM 
attenuation. The UV LF estimated with the method is 
fairly consistent with the bright end of the observed UV 
LF at z — 5.7, where the incompleteness is not severe. 
The bright end of the rest-UV LF down to 3a limiting 
magnitudes in the z' band can reasonably be reproduced 
from the observed Lya LF and the EWo-UV luminosity 
relation. However, the shape of the derived UV LF is 
closer to a power law than to a Schechter function. When 
we fit the UV LF to a Schechter function, the faint-end 
slope was found to be a = —2.4, which is much steeper 
than those derived from LBGs (Bouwens et al. 2007). 

The critical number density of ionizing photons neces- 
sary to keep the intergalactic hydrogen ionized was given 
in Maddau, Haardt, & Rees (1999). The corresponding 
critical star formation rate density (SFRD), p*(z), can 
be written: 

W.) = 0.013/,- (g) (l±i)' (H) 2 M 0yI -Mpc1S, 

where / esc is the escape fraction of ionizing photons, and 
C is the ionized hydrogen clumping factor of the IGM. 
Here, the authors assumed a Salpeter IMF, solar metal- 
licity, and SED at < 912A based on the Bruzual & Char- 
lot population synthesis model. Both / csc and C30 are 
highly uncertain. The / esc was generally estimated to 
be < 0.1 (Leitherer et al. 1995, Deharveng et al. 2001, 
Fernandez-Soto et al. 2003, Malkan et al. 2003, Shap- 
ley et al. 2006) for various galaxies at various redshifts. 
Iwata et al. (2009) and Inoue et al. (2011) obtained much 
higher / osc of <~ 0.5 for some LAEs at z = 3.1, and 
Vanzella et al. (2010) found one LBG at z = 3.795 with 
a direct detection of the Lyman continuum flux, which 
gave /esc > 0.15. Bouwens et al. (2010a) recently argued 
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for the possibility of high / csc at z ~ 7 to account for 
the observed steep UV- continuum slope without a nebu- 
lar emission contribution (see also Taniguchi et al. 2010). 
We assumed that the escape fraction did not depend ei- 
ther on luminosity or on redshift; though some models 
predict that it could increase up to ~ 0.8 in lower-mass 
galaxies at higher z (Ricotti & Shull 2000, Wise & Cen 
2009, Razoumov & Sommer-Larsen 2010). Siana et al. 
(2010), based on deep HST far-UV images, suggest that 
the escape fraction significantly increases from z = 1.3 
to 3. The IGM clumping factor was estimated to be 
C ~ 30 at z = 5 based on a numerical simulation, though 
it could be lower if higher density regions are less ionized 
(Gnedin & Ostriker 1997). A more recent simulation by 
Pawlik et al. (2009) suggests a much lower clumping fac- 
tor, C = 6. Bolton & Haehnelt (2007) suggest C < 3 at 
z — 6, based on measurements of the metagalactic pho- 
toionization rate, combined with a model for the ionizing 
photon mean free path. Figure 14 presents a comparison 
between for different parameter sets of (/ csc , C) in 
Equation (2) and the SFRD derived by integrating the 
UV LF estimated in this section. Assuming the complete 
escape of ionizing photons (/ osc = 1) and a homogeneous 
IGM (C = 1), the universe can be easily ionized only 
by bright LAEs at z — 6.5; however, in the more real- 
istic case with (/ esc < 1 and C » 1), the contribution 
of low-luminosity LAEs is important to the reionization 
process. 

Our method assumed the EWo-UV luminosity relation, 
which is found to be quite sensitive to the SFRD esti- 
mate. The dashed line above the solid line in Figure 14 
shows the SFRD estimates when w was reduced to half. 
The dashed line below the solid line in Figure 14 shows 
the SFRD estimates when the EWo = 15A cutoff was ap- 
plied instead of 10A. Both of these estimates are appar- 
ently inconsistent with the reliably observed bright end 
of the SFRD; therefore, only narrow ranges of parameters 
in the EWo-UV luminosity relation are acceptable. As 
mentioned in Section 5, an EWo larger than ~ 300 A is 
hard to explain using the usual Population II synthesis. 
When making a trial cutoff at the high-EWo end of the 
EWo-UV luminosity, the SFRD steepens sharply at the 
faint end. We also fit a one-sided Gaussian instead of an 
exponential to the EWo distribution function, but the 
estimate did not change significantly. The most uncer- 
tain parameter in the estimate was the faint-end slope of 
the Lya LF, a, which, thus far, was poorly constrained 
by observation. The magenta and cyan solid lines in Fig- 
ure 14 show the SFRD estimates when changing a from 
— 1.5 to —1.7 and —1.3, respectively. Although these two 
estimates hardly affect the bright end, maintaining con- 
sistency with the observations, they diverge strongly at 
the faint end. The steeper faint-end slope makes a larger 
contribution to the photon budget, as expected. 

The black solid line in Figure 14 shows the SFRD es- 
timate based on the UV LF of LBGs at z ~ 6 (Bouwens 
et al. 2007). Although the SFRD of LAEs seems to be 
larger than that of LBGs at the bright end, they are 
within the uncertainties (see the shaded region of Fig- 
ure 9 indicating the range of the UV LF of the LBG 
sample at z ~ 6 determined by various studies). At the 
faint end, the SFRD of LAEs increases more significantly 
than the SFRD of LBGs toward faint magnitudes, as ex- 
pected from the fact that the LAEs have steeper UV LF. 
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Fig. 14. — Comparison between p* for different parameter sets of 
(/csc, C) in Equation (2) and the SFRD (green solid line) derived 
by integrating the UV LF of LAEs at z = 5.7 estimated in Section 
7. The SFRD estimates based on the observed UV LF appear as 
the blue (red) shaded region for z = 5.7 (6.5), though these are 
incomplete at the faint end below the 3<r limiting magnitude. The 
dashed line above the solid line shows the SFRD estimates when ui 
is reduced to half. The dashed line below the solid line shows the 
SFRD estimates when the cutoff EWo = 15A instead of lOAis ap- 
plied. The magenta and cyan solid lines show the SFRD estimate 
when a is changed from —1.5 to —1.7 and —1.3, respectively. The 
black solid line shows the SFRD estimate based on the UV LF of 
LBGs at z ~ 6 (Bouwens et al. 2007). The minimum SFR inferred 
by the simulation by Nagamine et al. (2010) and approximate de- 
tection limit of JWST are shown by dotted lines. 

It should be noted that this result strongly depends on 
the faint-end slope, a, of the Lya LF at z — 5.7, which is 
poorly constrained by this study. Here, we provide only 
estimates for values of a = —1.3, —1.5 and —1.7. At low- 
z, Gronwall et al. (2007) derived a faint-end slope of Lya 
LF of a = -1.49l°;g, and Cassata et al. (2011) found 

a = -1.60tq ll at z~ 2.5 and a = -1.78±g;^ at z ~ 4. 
It is interesting that a steeper faint-end slope of the UV 
LF than of the LBG is predicted in Figure 14 even when 
a = — 1.3, which is shallower than those at low z. Stark 
et al. (2011) also suggested a higher faction of LAEs in 
their LBG sample towards fainter UV magnitudes. Con- 
sequently, the relative contribution to the photon budget 
of the LAEs compared with the LBGs increases toward 
fainter magnitudes. It finally exceeds the contribution 
of the LBGs at the faintest magnitude, which seems un- 
likely. Nagamine et al. (2010) employed a cosmological 
SPH simulation to reproduce some observational prop- 
erties of LAEs and found a critical threshold in stellar 
mass around M* ~ 10 7 Mq, below which star formation 
rapidly drops, as might be expected from the Kennicutt 
law. Therefore, we assumed that the minimum SFR 
that galaxies should have is 10 _2 M Q yr _1 , which cor- 
responds to the critical stellar mass in the simulation by 
Nagamine et al. (2010). When assuming a — —1.5, the 
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LAE contribution exceeds the LBG around Mjjv ~ —15. 
This might suggest that some of our assumptions are 
not reasonable; otherwise, the faint-end slope of Lya LF 
of LAEs should be shallower than —1.5. For example, 
given a = —1.3, the LAE contribution does not over- 
whelm the LBG contribution down to SFR= 1O _2 M 
yr _1 . However, we note here again that the derived UV 
LF of LAEs is based on many assumptions. For exam- 
ple, the faint-end slope of the SFRD becomes flatter even 
with a = —1.5 if much larger EWo are allowed at faint 
UV luminosities. On the other hand, we note that the 
LBG estimate is also based on an extrapolation of the 
observed UV LF to very faint levels. 

Despite many uncertainties, this might be the first es- 
timate of the contribution of the LAE population to the 
reionizing photons. The SFR estimate based on Lya lu- 
minosity is much more strongly affected by extinctions 
of dust, ISM, and IGM than the SFR estimate based 
on rest-UV luminosity. Our result indicates that low- 
luminosity LAEs could contribute significantly to the 
photon budget necessary for reionization, though these 
faint LAEs are below the current detection limits. Insofar 
as a universal EWo-UV luminosity relation is assumed, 
a steep faint-end slope in the UV LF is predicted even 
when the faint-end slope of the Lya LF is shallower than 
those at low z. The faint end of the LF of these ionizing 
sources is critical to the conclusion, and an accurate de- 
termination of the faint end slope of Lya LF will enable 
a precise estimate of the total ionizing photon density 
emitted by LAEs at this epoch. Bouwens et al. (2007) 
concluded that the number of i dropouts at z ~ 6 ap- 
pears to be approximately consistent with the numbers 
necessary to reionize the universe. Recent studies using 
HST/WFC3 identified more and fainter galaxies at even 
higher z, suggesting a faint-end slope that is sufficiently 
steep (Oesch et al. 2010, Bouwens et al. 2010b) to fully 
reionize the universe. The LAE sample based on NB 
searches exploring only a small redshift coverage is com- 
plementary to the estimate because LAEs should consti- 
tute a very young population at a single epoch among 
heterogeneous high-z star-forming galaxies sampled by 
the i-dropout method. Moreover, our LAE sample has 
the advantage of many spectroscopic confirmations, pro- 
viding an accurate estimate of sample completeness and 
contamination, which also increases the reliability of our 
photometric LAE sample. 

8. SUMMARY AND DISCUSSION 

We obtained extended spectroscopic confirmations of 
LAEs at z = 6.5 and 5.7 in the SDF, and our conclusions 
can be summarized as follows: 

1. We provided new identifications of 28 and 20 LAEs 
at z = 6.5 and 5.7, respectively. The discrimination of an 
LAE from other nearby emitters was based on a quanti- 
tative line asymmetric estimator, weighted skewness, S w . 
The total number of spectroscopically confirmed LAEs in 
the SDF is now 45 (54) at z = 6.5 (5.7), which means 
that 90% (74%) of the photometric candidates have been 
followed by spectroscopy. Our long campaign of follow- 
up spectroscopy shows that our photometric LAE sample 
is highly reliable, with low incompleteness and little con- 
tamination. 

2. We made more careful measurements of Lya lumi- 
nosity, both photometrically and spectroscopically, than 



in our previous study to determine Lya and rest-UV LFs 
more accurately. The non-square transmission curve of 
filters was taken into account in the photometric mea- 
surements, and we corrected for slit loss in the spectro- 
scopic measurements. The Lya fluxes measured photo- 
metrically and spectroscopically agree very well, showing 
that our measurement of Lya luminosity is accurate. 

3. We derived Lya LFs of LAEs at z = 6.5 and 5.7. 
With a large number of spectroscopic confirmations of 
our LAE sample and more careful measurements of lumi- 
nosities, the Lya LFs at both redshifts are more sharply 
determined than those in our previous studies. The sub- 
stantially improved evaluation of Lya LF of z = 6.5 
shows an apparent deficit from z = 5.7, at least at the 
bright end and a possible decline even at the faint end, 
though small uncertainties remain. 

4. We derived rest-UV LFs of LAEs at z = 6.5 and 5.7. 
We confirmed the result of our previous study, namely 
that the rest-UV LFs at z — 6.5 and 5.7 agree with each 
other, which is in clear contrast to the difference seen in 
the Lya LF, though the measurements of UV LF of LAEs 
still have large uncertainties. The rest-UV LF of LAEs 
is almost consistent with that of LBGs within the errors, 
at least at the bright end, suggesting that the fraction of 
LAEs among LBGs at z ~ 6 is higher than that at z = 3. 

5. The EWo distribution of our z = 5.7 LAE sample 
is almost consistent with other low-z samples, whereas 
the EWo at z = 6.5 seems to be systematically smaller 
than z = 5.7. This may be understood from the resultant 
LFs of these two epochs: UV-continuum LFs are almost 
identical, whereas the Lya LF at z — 6.5 is deficient 
compared to that at z = 5.7. There remains, however, 
the possibility that the trend is attributed to the artifacts 
of using BB band including the Lya line at z — 6.5. 
The EWo distribution at z = 6.5 shows an extended tail 
toward larger EWo compared with z = 5.7, though EWo 
measurements are not reliable when not detected in the 
BB band. The EWo and rest-UV luminosity relation 
shows an apparent deficit of LAEs with large EWo at 
bright UV magnitudes, and the upper-limit of the EWo 
increases with lower UV luminosity. 

6. We found no significant difference in the composite 
Lya line profile between z = 6.5 and 5.7, though the pos- 
sible velocity offset of Lya line from the systemic redshift 
might dilute the composite profiles. 

7. We tried to recover the rest-UV LF of LAEs at the 
faint end, where incompleteness is severe, by assuming a 
universal EWo-UV luminosity relation. When choosing 
reasonable parameters to fit to the observed EWo-UV 
luminosity relation, the bright end of the rest-UV LF 
at z = 5.7 was well reproduced from the observed Lya 
LF. Integrating this rest-UV LF permitted the first ex- 
perimental estimate of the photon budget of LAEs for 
reionization. The derived UV LF suggested that the 
fractional LAE contribution to the photon budget among 
LBGs significantly increases toward fainter magnitudes. 
Low-luminosity LAEs could contribute significantly to 
the photon budget, though this depends on the poorly 
constrained faint-end slope of the Lya LF. 

8.1. Implications for Reionization 

We confirmed our previous result that the Lya LF at 
z = 6.5 declines from those at z < 5.7 based on our 
deep LAE samples with a high spectroscopic identifica- 
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tion rate. This result is, more or less, consistent with 
Ouchi et al. (2010) and Hu ct al. (2010). The decline in 
the Lya LF of LAEs from z = 5.7 to 6.5 could be caused 
by the evolution of some intrinsic property of LAEs or 
by evolution in the ionization state of the IGM. It will be 
almost impossible to distinguish between these two pos- 
sibilities from the Lya LF alone; however, the rest-UV 
LF, which is not sensitive to the neutral IGM, may pro- 
vide an additional diagnostic of the cause of the observed 
trend. Interestingly, our measurements of the rest-UV 
LFs agree well within uncertainties between these two 
epochs at least at the bright end, where an apparent 
difference is observed in the Lya LF. This result might 
support the interpretation of neutral IGM attenuation. 
Assuming a fully ionized IGM at z = 5.7, the observed 
difference in the Lya LF suggests xui ~ 0.38 at z = 6.5 
based on the model of McQuinn et al. (2007). Another 
important aspect of ascribing the difference in Lya LF 
to the galaxy evolution of LAEs is that it should simul- 
taneously explain the difference between z = 5.7 and 6.5 
and the lack of evolution in the Lya LF from z = 3 to 
6. The model of Kobayashi ct al. (2010), which takes 
into account a possible evolution of / esc and reasonably 
reproduces Lya LF, UV LF, and the EW distribution of 
LAEs from z = 3.1 to 6.5, favors Xm ~ 0.4 at z = 6.5. In 
the model of Kobayashi et al. (2010), the observed large 
decline in the Lya LF at z — 6.5 could be partly caused 
by galaxy evolution; however, IGM attenuation is still 
required to account for the entire decline. The predicted 
value of the neutral fraction at the moment is strongly 
model-dependent. Dayal et al. (2008) concluded that the 
Lya LF at z = 6.5 could be reproduced by their model 
even with small xjji = 3 x 10~ 4 , though a strong increase 
in dust content was required to match the model to the 
Lya LF at z < 5. The amplitude difference in the Lya 
LF is still relatively small between these two epochs, sug- 
gesting that the universe was still largely (> 60%) ion- 
ized at z = 6.5. Such an amplitude difference is smaller 
than the uncertainties of the rest-UV LFs, though we 
concluded that the rest-UV LFs between the two epochs 
agree well at the bright end. More accurate measurement 
of the rest-UV LFs based on the larger LAE sample in 
the future will provide a stronger constraint on the Lya 
test. 

8.2. Cosmic Variance Uncertainties 

The decline in Lya LF from z = 5.7 to 6.5 might be 
caused by cosmic variance. The sample size is not yet 
large enough to statistically refute this interpretation; 
however, we found no difference in the rest-UV contin- 
uum LFs, which should be affected by cosmic variance 
from z = 5.7 to 6.5, if it exists. Ouchi et al. (2010) re- 
cently carried out an LAE survey at z = 6.5 with an 
FOV five times wider than this study. Though their 
number of spectroscopic confirmations is small, the re- 
sult is statistically robust and less affected by cosmic 
variance than ours. They concluded that the Lya LF at 
z = 6.5 shows a decline of 30% in Lya luminosity from 
z = 5.7, whereas the present study finds a 24% decline. 
Interestingly, they detected an apparent difference in the 
LF even at the faint end. Their LF difference alone con- 
strains the neutral fraction to xhi < 0.2^Q 2, though they 
did not detect the significant enhancement of clustering 
amplitude expected during reionization. Hu et al. (2010) 



also completed their LAE surveys at z = 6.5 and 5.7, us- 
ing an FOV five times wider than this study. Though the 
shape of their LF differs from this study and Ouchi et al. 
(2010), they also found that the Lya LFs differed based 
on their own estimates at these two epochs, indicating 
a 44% decline in 0*, whereas this study found a 34% 
decline. It is interesting that both of these wide-field 
studies confirmed a difference in the Lya LF between 
z = 6.5 and 5.7. Nakamura et al. (2011) also derived the 
Lya LF at z — 6.5, with a low number density only x0.3 
of our previous study. The degree of difference in the 
LF found in the present study, Ouchi et al. (2010), and 
Hu et al. (2010) are slightly inconsistent, possibly sug- 
gesting ficld-to-ficld cosmic variance. Otherwise, such a 
variance could be caused by a patchy reionization pro- 
cess, which has been suggested in the GP trough mea- 
surements showing a substantial variation in IGM trans- 
missions among different QSO lines of sight around z ~ 6 
(Djorgovski et al. 2006). Such a process is expected to 
be caused by primordial clustering or an initial large- 
scale structure of ionizing sources. Much wider surveys 
of high- z LAEs arc required to obtain statistically signif- 
icant constraints on any variance. 



8.3. Photon budget of Reionization 

In Section 7, we attempted to estimate the rest-UV LF 
of LAEs using the Lya LF, assuming a universal EWo- 
UV luminosity relation. Interestingly, this estimate re- 
produces the bright end of the observed UV LF reason- 
ably well, enabling an estimate of the LAE contribution 
to the photon budget required for reionization. We note, 
however, that this method has large uncertainties, in- 
cluding the universality of the EWo distribution function 
and the EWo-UV luminosity relation. Pentericci et al. 
(2010) indicated that a correlation between Lya strength 
and age or SFR might change with cosmic time, poten- 
tially changing the EWo distribution function. Nilsson 
et al. (2009) suggested no correlation between EW and 
UV luminosity; however, this assumption, coupled with 
the observed Lya LF, produces a much steeper faint-end 
slope in the UV LF of LAEs, exacerbating the incon- 
sistency with the UV LF of LBG. The universality of 
the EWo distribution and EWo-UV luminosity relation 
should be verified at lower z based on a larger LAE sam- 
ple. The observed large scatter in Lya EWo at faint 
UV luminosity could be due to small amounts of dust 
extinction (Verhamme et al. 2008) and possibly to a 
stochasticity of their complicated duty circle (Nagamine 
et al. 2010), which has not yet been clearly revealed. Fur- 
thermore, more accurate determinations of the faint-end 
slope of the Lya LF, / osc , and C are required to pro- 
vide a stronger constraint on the photon budget. The 
James Webb Space Telescope (JWST) will extend the 
current observational limit down to htlab ~ 31, which 
corresponds to Mjjy = —15.6 at z = 5.7 (see Figure 14). 
Significant detections of LAEs down to this faint end will 
allow a more precise estimate of the photon budget. The 
JWST/NIRCam and TFI will catch LAEs at higher z 
far beyond the current frontier of the distant universe. 
The derived Lya LF during the early reionization phase 
will be more sensitive to the neutral fraction, giving a 
stronger constraint on the history of reionization. 
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